Fractional pumping of energy into a ballistic quantum ring. 
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We consider the energy stored in a one-dimensional ballis- 
tic ring with a barrier subject to a linearly time-dependent 
magnetic flux. An exact analytic solution for the quantum 
dynamics of electrons in the ring is found for the case when 
the electro-motive force £ is much smaller than the level spac- 
ing, A. Electron states exponentially localized in energy are 
found for irrational values of the ratio A = A/2e£. Relax- 
ation limits the dynamic evolution and localization does not 
develop if A is sufficiently close to a rational number. As a 
result the accumulated energy becomes a regular function of 
A containing a set of sharp peaks at rational values (fractional 
pumping) . The shape of the peaks and the distances between 
them are governed by the interplay between the strength of 
backscattering and the relaxation rate. 
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Physical properties of mesoscopic systems are governed 
by quantum interference. Several phenonema of such a 
nature have been discussed for systems close to equilib- 
rium. Persistent currents in multiply-connected systems 
|ij as well as universal fluctuations of the conductance 
are important examples j^j . Coherent dynamics remains 
crucially important in situations far from equilibrium 
provided the energy associated with the phase breaking 
rate is less than the characteristic rate of redistribut- 
ing electrons in energy space. Consequently, one can ex- 
pect pronounced mesoscopic behavior even in strongly 
biased mesoscopic devices, where the dynamics can be 
effectively tuned by external electric or magnetic fields. 

In this paper, we consider an example of such a sys- 
tem, namely a single-channel mesoscopic ring subjected 
to a nonstationary perpendicular magnetic field, linearly 
dependent on time. We concentrate on the energy ac- 
cumulation in such a system. To investigate the role of 
interference, we take into account electron backscatter- 
ing from a single potential barrier, embedded in the ring. 
Tuning the transmission through the barrier by gate po- 
tentials one can influence the interference pattern and in 
this way significantly change the dynamics. 

Impure onducting rings have been extensively dis- 
cussed in connection with energy dissipation in meso- 
scopic metallic systems. Gefen and Thouless |3|,|| have 
suggested that disorder leads to a localization of elec- 
trons in energy space. Consequently, a time-dependent 
magnetic flux, 4>, through the ring could induce a (dc) 



circular current only in the presence of phase breaking 
processes. This work was continued in Ref. Q by a nu- 
merical analysis of the role of dissipation. A numerical 
analysis of localization in a dissipation-free impure ring 
was performed in Ref. ||. 

The vanishing of the flux-induced current in an im- 
pure dissipation-free ring ||,|| is drastically different from 
what happens in a perfect ballistic ring, where a non- 
equilibrium flux driven current diverges when the dissi- 
pation goes to zero. This makes the crucial importance of 
the backscattering strength clear: By tuning the height 
of the barrier one can cross over from one regime to the 
other, and in this way control the pumping of energy into 
the system. This issue, not addressed in previous work, 
is the subject of the present paper. 

We will show that the scenario of the cross over is 
as follows. Consider the conductance of the ring, G, de- 
fined as the ratio between the circulating current and the 
electro-motive force £ = — (l/c)$ induced in a ring of ra- 
dius ro by a magnetic field linearly dependent on time. If 
scattering is strong, G oc £~ 2 . As the scattering strength 
decreases, a set of peaks in the function G — G{£) ap- 
pears. The peaks correspond to rational values p/q of 
the quantity A = A/2e£, where A = h 2 Np/mrQ. Here 
Np is the number of filled electron states while m is the 
effective mass. The shape of the peaks and the distances 
between them are governed by an interplay between the 
height V of the potential barrier and the relaxation rate, 
v, the maximum value of q being determined by the con- 
dition r q /q ~ hv/e£. Here r = exp(— £ c /£) is the ef- 
fective amplitude of Zener tunneling through the energy 
gaps in the electron spectrum, £ c = n 2 V 2 /2Ae. The 
peak structure near maximum can be described by the 
interpolation formula 
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Here g = G/Gq, Go — e 2 /h, while T)(e) is a smooth func- 
tion of e. If |e| < r q /q 2 the function r\ ~ 1, beyond this 
region it decreases as |e| increases. As the barrier be- 
comes more transparent, r — ► 1, the inter-peak distance 
(determined by the maximum value of q) decreases. Fi- 
nally, the peaks overlap forming the ballistic-like conduc- 
tance g = A/hv. 

To understand the result conjectured above let us con- 
sider the electron energy levels in the vicinity of the Fermi 
level, Ef- Here the energy dispersion can be consid- 
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ered to be linear. In a ballistic ring, one has then two 
sets of adiabatic energies Ei(&) corresponding to clock- 
and counterclockwise motion (Fig. 1). The sets coin- 
cide and scattering from the barrier opens gaps for those 
flux values where E\ — Ep + /A/2, I = 0, ±1,±2, . . .. 
Consequently, the energy pumping into the system by a 
slowly varying magnetic flux can be mapped to the one- 
dimensional motion of a quantum particle in the field of 
periodically placed scatterers (cf. Refs. Landau- 
Zcncr tunneling (with the amplitude r introduced above) 
through the gaps corresponds to forward scattering while 
reflection from the gaps is similar to backscattering. The 
important difference from the usual impurity problem 
is that there is no translational invariance at an arbi- 
trary value of the driving force £. This invariance is 
only present for rational values p/q of the dimension- 
less quantity A In this case we arrive at a super- 
lattice containing q "impurities" per unit cell. As a re- 
sult, the motion along the E-&xis is described by q al- 
lowed bands, the "velocity" being ve = E ~ At 9 /to 
(here to = h/e£ is twice the time interval between two 
sequential Landau-Zener scattering events). Since the 
upper bound of the Brillouin zone is Airh/qA, the cor- 
responding bandwidth for motion along the _E-axis is 
W ~ VE^ith/q/S. ~ A-KTiT q /t q. At rational values p/q 
of the quantity A the electron rotates around the ring 
an integer number, 2p, of times during the time the en- 
closed magnetic flux changes by q flux quanta. As a re- 
sult motion along the _E-axis can be mapped onto motion 
in a one-dimensional periodic potential, the correspond- 
ing eigenstates being extended. If p/q is irrational the 
motion along the iS-axis is equivalent to motion in a one- 
dimensional quasi-periodic potential. It turns out that 
the corresponding states are then localized (see below) 
in spite of the fact that there is no real disorder in the 
system. The localization radius, Ri oc , can be estimated 
for A = p/q + e, \e\ <C 1/q as follows: In this case the 
phase mismatch due to finite e can be thought of as be- 
ing due to a quasicassical potential U(E) = eaE with 
a = &itTi/ Ato- One can see that a finite e gives rise to 
band bending which creates semiclassical turning points 
for the propagating modes along the E-&xis. The local- 
ization length is in fact one half of the distance between 
the turning points produced by the upper and lower band 
edges, Rioc ~ W/2a\e\ = Ar q /Aq\e\. Consequently, the 
localization time is ii oc ~ ^Rioc/ve ~ W^M- 

The manifestation of localization in the energy pump- 
ing depends on the product vt\ oc . When vt\ oc 3> 1 local- 
ization has not developed and the band picture of energy 
pumping is relevant. The conductance is estimated as (cf. 
Rcf. |J) G = P/£ 2 , where P is the average energy accu- 
mulation rate. The quantity P, in its turn, is determined 
as v8E ■ N(5E). Here 5E ~ ve/v is the energy accumu- 
lated by a single state, while N(SE) ~ 8E/A is the num- 
ber of involved states. It follows that g ~ T 2q (A/Tiv). 
If vt\ oc « 1, on the other hand, G is determined by 



hops between intra-band localized states. In this case, 
SE ~ 2R loc , and we obtain g ~ TivAT 2q /(e£qe) 2 . These 
estimates are consistent with the first term in Eq. ([j]). 

The consideration above are valid if the localization 
length is greater than the unit cell size, 2Ri oc ^> qA/2. 
In the opposite limit the band picture fails and the veloc- 
ity is dominated by inter-band transitions. In this case, 
SE ~ qA, N(SE) ~ q, and g ~ q 2 TivA/{e£) 2 . This 
estimate corresponds to the second term in Eq. ([!]). 

The following model is employed. The electron system 
is described by the Hamiltonian 

f d t\ ch 
H(t, <p) = -Aa z [i — + - + Wtopfa), to = — r (2) 
V d<P W |e|$ 

Here <7i are Pauli matrices. We are interested in the 
current, averaged over the time to, 

i /-t+to/2 

J(t) = - / dt'Trpj j=eip. (3) 

tO Jt-to/2 

The single-electron density matrix, p, is calculated from 
the equation 

^ = jr\p,H]-v(p-fo(H{t))), (4) 

where fo is the Fermi function. The formal solution of 
Eq. (||) can be expressed in terms of the evolution opera- 
tor u(t, t') for the systems described by the Hamiltonian 




p(t) = is[ t dtV^'-^OMwCO)^ (*,*')• (s) 

J — oo 

We are interested in the case of weak scattering, i.e. when 
the matrix element V of the scattering potential cor- 
responding to a momentum transfer of ~ 2pp is much 
smaller than the inter- level spacing A. In such a sit- 
uation, the impurity potential is important only in the 
vicinity of the times when the adiabatic levels for clock- 
wise and counterclockwise motion cross. It creates gaps 
near the crossing points. Consequently, one can discrim- 
inate between intervals of ballistic evolution (duration 
~ to/2) and intervals of Landau-Zener scattering. The 
typical duration of such an interval is < WtoH/A (cf. 
Ref. ||). Consequently, at Ato/K ^> 1 the Landau-Zener 
scattering is strongly confined within narrow intervals 
and can be described in terms of the scattering matrix 

S = cxp(i6») [pexp(i^cr 2 ) + ira x ] , p = y/l-r 2 (6) 

where t has been introduced above. It turns out that 
the physical values of interest here do not depend on the 
phases 9 and £. For simplucity we put them equal to 
zero. The expression for J has the form 

Y = -t 2 + J2 e~ 2£im Tr f- m JT m T , (7) 
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where Jq — eA/h is the amplitude of the persistent cur- 
rent, v = vto/2, T = exp(— itp)u(to, 0), 

I /-to/2 

J = — dtu + (t,0)a z u(t,0) , 

*o J -t /2 



^««W M * (,)) «M) 

-to/2 dt 



Here we have employed the relation u(t + mto, t' + mto) = 
exp(im(p)u(t,t') exp(—imtp) which immediately follows 
from a similar symmetry property of Ti. 

The problem of calculating the current is reduced to an 
analysis of the unitary operator T. Having in mind the 
periodicity in tp of all quantities we introduce the basis 



e ±i(N F +n) V / ^ 

' n ' s = ± ^ ' s± ' s+= (o 



s_ = 



In this representation, the operator T can be expressed as 
a direct product of operators acting in n and pseudo-spin 
(s) spaces, 



T = pSe l ™ a * 



•R 



±e 



where the operators R± and n are defined as R±\n, s) — 
|n=Fl, s), and h\n, s) = n\n, s), and a is the fractional part 
of the quantity A introduced above, o~± = [cr x ± ia y )/2. 
In the n, s-representation the operators J and T have the 
form J = V5 ntn >, f = V5 0tn 5 ,n>, V = r 2 a z - pra y f. 

It is straightforward to show that the unitary operator 
T has the following properties: 

R-TR+ = e - il7ra f, a y f*(-h)a y = f(n) 

These result in the following relations between the 
eigenstates ipp{n) of the operator T having eigenvalues 
exp(i/3): 



ipp(n + 1) 
a y tpp(-n) 



(8) 



Hence, the spectrum of T can be expressed in the form 
exp(z/3^) where = ±/3 (a) — 4irar, r = 0,±1,±2, . . . 
One can prove that the properties (||) allow one to gener- 
ate a complete set of eigenstates provided i/)p is known. 
Furthermore, the vector equation for the eigenstates of 
T reduces to a scalar equation 
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B(m) = 0, 



where B is a linear combination of the components ipg. 
Its solution allows one to determine both the spectrum 
and the wave functions. The detailed analysis will be 
published elsewhere ||. 

The results are different for the cases of rational and 
irrational values of a. 



a. Irrational a. The analysis shows that (3q = it a and 
the set of {r} is infinite. The eigenfunctions have the 
form 



„^inan(2n±l) f2n / i x (k) 

^ a = - = / dke- 2ikn '' 



2tt 



i) 



ix(k—Tza)+ik J ' 



x(*) = -E 



i=i 



t 1 cosl(k-w/2) 
I sin ixal 



(9) 



It can be shown that the infinite series (||) converges for 
almost all irrational values of a, and x(k) is an analytic 
function ]Tc| ] . Consequently, the eigenfunctions are expo- 
nentially localized, the localization length i?i oc in energy 
space being 



R loc 

A 2 



E 



4 sin nal 



Here |/3) = xpp(n) is the eigenstate of the operator T. 
One can see that in the vicinity of rational values p/q of a 
the localization length i?i oc diverges as Ar q /2irq\a — p/q\ 
in agreement with the qualitative estimates given above. 
The expression for the dimensionless conductance g at 
isto <C 1 yields 



g = nhv 



A (2R l 



{e£f 



b. Rational a{= p/q). Since the problem is transla- 
tionally invariant in n-space. Consequently, the eigen- 
states can be labeled by a quasi momentum /C. The 
spectrum is now given by 



a lv\ n P 2 • / ? • 1*" 
Po(/C) = arcsm < T q sm — 

q q { 2 



While the current can be expressed as expression 

- = + £ I dicing 



(10) 



±.r=0 



r v 2 + sin 2 $r 



*r = 2 ' ^rK^ 2^W(3oi Vc Ppt) ■ I 11 ) 

z=o 

Here <Ag± (I) is the Bloch amplitude corresponding to the 
eigenstate V'wcO- A n exact expression for the eigenfunc- 
tion ij)R Q (n) shows the limiting transition to the expres- 
sion (0) as q,p — * oo, p/q = const. The current calcu- 
lated according from Eq. ( pi] ) remains continuous also. 
Thus, Eq. ( [Ti"l ) with large enough q and p can be used as 
a good approximation for irrational a-values. The results 
of such a calculation is shown in Fig. 2. 

The following two assumptions have been implicitly 
made in our discussion: (i) the electron dynamics is gov- 
erned by a linear dispersion law; (ii) the energy gaps as 
well as the scattering matrix S do not depend on the 
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particular energy levels involved. Assumption (i) is valid 
if the number of involved states (limited by the relax- 
ation rate) is much less than Np. Estimates show that 
(i) is valid if (hv / A) 2 {A/ e£)N F > 1. The first factor in 
this product is small while the two othes are large; the 
criterion can be met under realistic experimental condi- 
tions. Assumption (ii) is valid if the Fourier component 
of the barrier potential, J V(<p)e 2mtp dp, is essentially n- 
independent for relevant n. This is the case if the scatter- 
ing potential is confined to a region of width 5p <C 1ixvt§. 
Note that the inequality vt$ <C 1 is essential for main- 
taining a noticeable energy pumping. 

Another approximation is that we have allowed for re- 
laxation in the simplest possible way by using a single re- 
laxation time [cf. Eq. (|)]. This is adequate if relaxation 
is caused by real space transfer of electrons between the 
ring and a surrounding reservoir. If the electron energy 
level spectrum in the reservoir is continuous the lifetime 
of an electron state in the ring with respect to this mech- 
anism is almost independent of quantum number. This 
mechanism allows us to describe electron states in the 
ring as pure quantum states with a relaxation time given 
by the time of decay through escape to the reservoir. 
The exact results obtained above are relevant for the case 
when such an "escape" mechanism dominates. Internal 
inelastic relaxation processes in the ring can lead to a sig- 
nificant difference between phase- and energy relaxation 
rates and requires a separate treatment. However, in the 
most interesting case of efficient Landau-Zener tunneling, 
the intrinsic inelastic processes must involve large mo- 
mentum transfers and are therefore strongly suppressed 

& 

Finally, it would be interesting to discuss the cross-over 
from the above picture to a case where more than one 
scatterer introduces disorder to the problem || . Interfer- 
ence between reflections from different scatterers would 
induce an energy dependence of the gaps in the spectrum 
and brake translational invariance in time. This would 
strongly affect the quasi-ballistic regime and presumably 
suppress the absorption peaks in Fig. 2. 

In conclusion the quantum electron dynamics problem 
in a single-channel ballistic ring with a barrier subjected 
to a linearly time-dependent magnetic flux has ben solved 
exactly. Exponential localization in energy space has 
been proven. Finally, we have shown that the dc-current 
has a set of peaks with fractional structure when plotted 
as a function of the induced electro-motive force. This 
structure is strongly sensitive to the barrier height, as 
well as to the relaxation rate. 
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FIG. 1. Diagram showing coincidence of flux-driven en- 
ergy levels (corresponding to clockwise and anti-clockwise mo- 
tion of electrons around ring) at a special flux value (cf. text). 




FIG. 2. The normalized current J/ Jo as a function of 
a = p/q for different Landau-Zener tunneling amplitudes, r. 
Jo = eA/h. A of v = 0.2 for the dimensionless relaxation rate 
was used. 
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